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Abstract 

Gliders in one-dimensional cellular automata are compact groups of non- 
quiescent and non-ether patterns (ether represents a periodic background) 
translating along automaton lattice. They are cellular- automaton anal- 
ogous of localizations or quasi-local collective excitations travelling in a 
spatially extended non-linear medium. They can be considered as binary 
strings or symbols travelling along a one-dimensional ring, interacting with 
each other and changing their states, or symbolic values, as a result of in- 
teractions. We analyse what types of interaction occur between gliders 
travelling on a cellular automaton 'cyclotron' and build a catalog of the 
most common reactions. We demonstrate that coUisions between gliders 
emulate the basic types of interaction that occur between localizations in 
non-linear media: fusion, elastic coUision, and soliton-like collision. Com- 
putational outcomes of a swarm of gliders circling on a one-dimensional 
torus are analysed via implementation of cyclic tag systems. 

Keywords: cellular automata, particles, travelling localizations, coUi- 
sions, beam routing, universality 



1 Introduction 

The era of unconventional computers — implementations of computing schemes 
in physical, chemical and biological substrates, and interpretation of the be- 
haviour of naturai systems' in terms of computation - has brought us a plethora 
of originai prototypes of future and emerging computing architectures. Some ex- 
amples of this new paradigm are polymer- and conductive foam based extended 
analog computers [42^ , chemical reaction- diffusion computers [2TJ [4] , aromatic 



molecular computers [10] , slime mould computers [50] , micro-fiuidic based com- 
puters [20], enzymatic logicai circuits molecular arithmetical circuits [58] 
to name but a few. Many of these novel computing devices work on the princi- 
ple that Information is represented by localised states of naturai systems (phase 
or diffusive waves, propagating pseudopodia, electron density, conformation) 
which can then be represented by cellular automata or other discrete automata 
models (at least in principle) [1]. Most models of unconventional computers 
suffer, up to different degrees, from boundary problems. The systems are con- 
fìned to some experimental arena, e.g. test tube or a Petri dish, or even be 
controlied and programmed by externally imposed geometrical constraints (e.g. 
channels in excitable chemical computers). It would be incredible useful to pro- 
duce a model which is self-contained, is not dependent on external constraints 
and which can, subject to resources and energy supplied, function indefìnitely. 
A particle "super collider" may be the best candidate for such a universal model 
of boundary free computation. 

In his seminai paper "Symbol super colliders" Tommaso Toffoli envis- 
aged far fetching analogies between physical implementations of particle super 
colliders, lattice gas and one-dimensional cellular automata. He suggested that 
the concept of symbol super collider - where myriad of tokens run along in- 
tersecting rings and interact with each other to produce new tokens - can be 
used in designing novel types of nature-inspired computing devices [TTl [59] . In 
present paper we develop Toffoli's ideas further and provide computational im- 
plementations of particle 'accelerator' or test rigs, implemented in elementary 
cellular automata. 

The paper is structured as follows. One-dimensional cellular automata, ahis- 
toric and with memory, are introduced in Sect.|2j The concepts of supercollider 
and ballistic computing are presented in Sect.jSj Section [4| shows how essential 
elements of collision-based computers can be implemented via glider interac- 
tions in one-dimensional cellular automata, and computational capacities are 
presented in Sect. [5j Outcomes of the presented results and directions for fur- 
ther studies are discussed in Sect. [6j 

2 One-dimensional cellular automata 

A cellular automaton (CA) is a quintuple (S, (/p, /ì, cq, L) based on a one-dimensional 
lattice L of cells, where each celi i G takes a state from a finite alphabet 
S. A sequence 5 G of n cell-states represents a string or a global configura- 
tion c on E. We write a set of finite configurations as E^. Cells update their 
states via an evolution rule (/:? : ^ S, such that /i represents a celi neigh- 
bourhood that consists of a centrai celi and a number of neighbours connected 
locally. There are different neighbourhoods and if k = |E| then we have 
different evolution rules. 

An evolution diagram for a CA is represented by a sequence of configurations 
{ci} generated by the global mapping $ : ^ S"^, where a global relation is 
given by as $(c^) c^+^. cq is the initial configuration. Celi states of a 
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configuration é are updated simultaneously by the evolution rule as: 

Lp{x\_^, . . . ,x\, . . . ,x\,^^) ^ x'r^^ . (1) 

where i indicates celi position and r is the radius of neighbourhood ji. Thus, 
the elementary basic CA class represents a system of order (/c = 2, r = 1) 
(in Wolfram's nomenclature [6ll), the well-known ECA. To represent a specifìc 
evolution rule we write name of the rule in a decimai notation, e.g. (^rhq. 

Conventional cellular automata are memoryless: the new state of a celi de- 
pends on the neighbourhood configuration solely at the preceding time step of 
if>. CA with memory are an extension of ECA in such a way that every celi Xi 
is allowed to remember its states during some fixed period of its evolution. CA 
with memory have been proposed originally by Alonso-Sanz 0[6l[7]. 

We implement a memory function (j) as follow: 

<P{x'-\---,x'-\x'i)^Su (2) 

where r < t determines the degree of memory and each celi G S is a state 
function of the series of states of the celi Xi with memory up r. To execute the 
evolution we apply the originai rule sl_i^ 5^, sl_^i^ . . .) ^ xl^\ Thus m 

CA with memory, while the mapping cp remains unaltered, historic memory of 
ali past iterations is retained by featuring each celi as a summary of its past 
states from cj). We can say that cells canalize memory to the map cp [7]. 

Let US consider the memory function ^ as a majority memory (pmaj ~^ 5^, 
where in case of a tie given by Si = Sq from (j) we take the last vaine Xi. So 
(l)maj function represents the classic majority function for three values as 
follows: 

(a A ò) V (ò A c) V (c A a) 

that represents the cells . . . , x-) and defines a temporal ring before 

to get the next global configuration c. Of course, this evaluation can be for any 
number of values. In this way, a number of functional memories may be used 
such as: minority, parity, alpha, etc. (see [7 ). 

Evolution rules representation in ECA with memory as given in [30] |3T] is 
as follows: 

4^CARm:T (3) 

where CAR is a decimai notation of a particular ECA rule and m is a kind of 
memory used given with a specific vaine of r. Thus for example, the majority 
memory (maj) incorporated in ECA Rule 30 employing five step of a cell's 
history (r = 5) is denoted as (j)R30maj:5' The memory is a function of the CA 
itself, see schematic explanation in Fig. [l] 
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ECA with memory (historic) 



t-T 



tempora! ring storing memory 
{S^} 



Figure 1: Illustration of cell-state transition in ECA (left) and ECA with mem- 
ory (right). 



3 TofFoli's supercoUider 

In the late 1970s Fredkin and Toffoh developed a concept of a general-purpose 
computation based on bahistic interactions between quanta of information that 
are represented by abstract particles [60]. The Boolean states of logicai vari- 
ables are represented by balls or atoms, which preserve their identity when they 
collide with each other. They came up with the idea of a billiard-ball model 
of computation, with underpinning mechanics of elastically colliding balls and 
mirrors reflecting the balls' trajectories. Later Margolus developed a special 
class of CA which implement the billiard-ball model. Margolus' partitioned CA 
exhibited computational universality because they simulated Fredkin gate via 
collision of soft spheres ^36j . 



u V u V u V 




(a) (b) (c) (d) 



Figure 2: Basic schemes of ballistic collision between localizations representing 
logicai values of the Boolean variables u and v. 

The following basic functions with two input arguments u and v can be 
expressed via collision between two localizations: 
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1. f{u^v) = c, fusion (Fig. |2^) 

2. f{u^v)=u^v^ interaction and subsequent change of state (Fig. ^p) 

3. fi{u,v) {u^v) identity, solitonic collision (Fig. [2]3); 

4. fr{u^v) ^ {v^u) reflection, elastic collision (Fig. |2]i); 

To map Toffoli's supercollider [60] onto a one-dimensional CA we use the 
notion of an idealized particle p G S+ (without energy and potential energy). 
The particle p is represented by a binary string of celi states. 




(a) (b) (c) 



Figure 3: Representation of abstract particles in a one-dimensional CA ring. 

Figure [3] shows two typical scenarios where particles pf and ps travel in a CA 
cyclotron. The first scenario (Fig. [3|i) shows two particles travelling in opposite 
directions which then collide. Their collision site is shown by a dark circle in 
(Fig. [3^). The second scenario demonstrates a typical beam routing where a 
fast particle pf eventually catches up with a slow particle Ps at a collision site 
(Fig. [3]3). If the particles collide like solitons [24], then the faster particle pf 
simply overtakes the slower particle Ps and continues its motion (Fig. [s]^). 

Typically, we can fìnd ali types of particles manifest in CA gliders, including 
positive negative and neutral displacements [44], and also composite 
particles assembled from elementary localizations. Let us consider the case 
where a quiescent state is substituted by cells synchronized together as an ether 
(periodic background). This phenomenon is associated with ECA Rule 110 
^i^iioQ Its evolution space is dominated by a number of particles emerging in 
various different orders, some of which are really quite complex constructions. 
Consequently, the number of collisions between particles is increased. Each 
particle has a period, displacement, velocity, mass, volume, and phase [44l [45]j^ 
Figure [4] displays a typical collision between two particles in (^rhq. As a result 
of the collision one particle is split into three different particles (for details please 
see [37]). The pre-collision positions of particles determines the outcomes of the 
collision. 

^Rule 110 repository ^http : //uncomp . uwe . ac . uk/genaro/Rule llO . html 

■^A full description of particles in Hule IIU is available at |http: //uncomp. uwe . ac .uk/| 
jgenaro/rule 1 10/glidersRule 110 . html| 
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Figure 4: Example of a particle collision in (^rhq. Particle collides with 
particle Pq giving rise to three new particles — p^, p^^, and p\^ — that are 
generated as a result of the collision. 

To represent particles on a given beam routing scheme (see Fig.[3|, we do not 
consider ether confìguration in (fRuo because this does not affect on collisions. 
Figure [5] displays a one-dimensional confìguration where two particles collide 
repeatedly and interact as solitonic so that the identities of the particles are 
preserved in the collisions. A negative particle p'^ collides and overtakes a 
neutral particle p'^^ . Figure 5a presents a whole set of cells in state 1 (dark 
points) where the ether confìguration makes it impossible to distinguish the 
particles: and p'^ . However, we can apply a fìlter and thereby select particles 
from their background ether (Fig. |5]3)|^ Space-time confìgurations of a cellular 
automaton exhibiting a collision between particles p'^ and p^ are shown in 
Fig#. 

Filters selected in CA are a useful tool for understand "hidden" properties 
of CA. This tool was amply developed by Wuensche in the context of automatic 
classifìcation of C A . The fìlters were derived from mechanical computation 
techniques [22], pattern recognition [57 , and analysis of cell-state frequencies 
[63] . Thus, a fìlter is a sequence of cells that have a high frequency in the 
evolution space. Such d-dimensional string repeat periodically, coexisting with 
any complex structure without altering or disturbing the global dynamics. 

^Ring evolution was simulateci with DDLab available inlhttp: //www.ddlab.orgl 
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(a) (b) (c) 

Figure 5: Example of a soliton-type interaction between particles in one- 
dimensional cellular automaton (pmio: (a)-(b) two steps of beam routing, (c) ex- 
act confìguration at the time of collision. 

4 Ballistic collisions in cellular automata 

In this section we analyze examples of ballistic collisions between particles in 
ECA and ECA with memory. Let us start with ECA with memory governed 
by the rule (j)R22maj-A There are only two types of particles in this rule 

^0K22maj:4 = {d L 1 9 r] ' Thclr properties are easy to infer: Both particles have a 
volume of a perfect square of 11 x 11 cells, a mass of 35 cells, and they translate 
2 cells per 11 time steps (or iterations of CA evolution). The particle ql has a 
negative slope with a velocity of — ^ and the particle qr has a positive slope 
with a velocity of 

The rule 4>R22maj-A supports two types of ballistic collisions: fi{u^v) and 
fr(u^v). Figure shows how a number of soliton interactions can be synchro- 
nized as a identity collision fi{u^v)^ where both particles can cross their own 
trajectories. This result can be achieved by selecting a particular phase of each 
particle as encoded by its initial condition. Different initial conditions lead to 
different reactions [38 . For example, in Fig. [gJ) we can see particles undergo- 
ing elastic collisions similarly to a lattice gas model [60 . This is a reflection 
collisions fr{u^ v). 



particle 


velocity 




2/2 = 1 


^ 


-2/2 = -1 


9o 


0/4 = 


9e 


0/4 = 


gun 


0/32 = 



Table 1: Properties of particles in rule (PR54 [32- 
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Figure 6: Ballistic coUisions in CA (l)R22maj-A- (a) identity or soliton coUision, 
and (b) refiections. 
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Figure [t] presents solitonic reactions (identity ballistic collision) for rules 
^R54 (memoryless ECA) and (t)mmaj-A (ECA with memory). The rule (^7^54 has 
been studied in detail in [TTl [22l |32l In representing particles in (fR^/^ as 

^^R5A — { ^7 ^07 gun}, we follow Boccara's et al. notation [11 . Table [l] 
gives the relation between particles and velocityj^ In this way, to produce a 
required reaction we must code a particular initial condition, where a lÈ particle 
is ready to collide with a fe" particle and both particles collide with the same 
phase with a stationary particle Qe (Fig. [7^). 

Figure [7|^b) shows an ECA with memory which preserves identity ballistic 
collision starting with any random initial condition. Thus, (pRQmaj-A evolves 
with two particles Qcf)RQ^^j.4^ = l"^, V}- The particles' properties are easy to 
calculate. The ^ particle has a volume of 5 x 6 cells, a mass of 12 cells, and 
moves 2 cells in 5 generations (positive slope). While V particle has a volume 
of 5 X 3 cells, a mass of 7 cells, and moves 2 cells in 5 generations (negative 
slope). 

As the example above demonstrates, by using ECA and ECA with memory 
we can experimentally study a variety of ballistic collisions. 



5 Beam routings and computations 

In this section we will exploit beam routing to produce some complex construc- 
tions that are based on particle-collisions. An additional effort is necessary to 
code initial conditions for every particle and recognize the most suitable phase 
for each particle in order to produce the desired reactions. Also, we will show 
how the beam routing can be used in design of computing based-collisions. 

Let US consider a ECA with memory, rule 4>RZQmaj:d> [30]. We want to 
implement a simple substitution function addToHead working on two strings 
wi = Al, . . . , and W2 = ^i, . . . , Bm^ where n, m > 1. For example, if 
wi — AAA, W2 = BBB and = W1W2 then addToHead(|i(;2|) means produce 
'^3 = W2W1. To implement such a function in ^i?30maj:8 every quantum of data 
is represented by a particle. Particles gi and ^2 are coded in order to reproduce 
a soliton reaction. The codifìcation is not sophisticated. However, a systematic 
analysis of reactions is required. A periodic gap and one fìxed phase between 
particles is suffìcient to reproduce addToHead function for any string A^B^. 

Figure [8] shows the evolution of (j^Rsomaj-.s starting from an initial condition 
coded by particles representing the string AAAAAAAABBBBBBBB. Using 
function addToHead we produce the final string BBBBBBBBAAAAAAAA af- 
ter 5,000 generations. The first snapshot (Fig.|8^) shows the initial configuration 
and the first 391 generations, the middle snapshot (Fig. ^) demonstrates how 
string wi acts on string W2 while preserving the Information (soliton or identity 
reaction), and the third snapshot (Fig. [s]^) shows the final global configuration 
(i.e. string W2W1 processed in parallel by (pRsomaj-.s)- 



Rule 54 repository in http://micomp.uwe.ac.uk/genaro/Rule54.html 

^Velocity is calculated as the number ot cells a particle displaces in one unit of discrete 
time. 
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(c) 

Figure 8: A simple substitution system processing the word A^B^ to B^A^ 
with (/>i?30maj:8- The final required state is reached at the 5,000th generation by 
synchronization of multiple soliton reactions. 
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Figure 9: Beam routing performing identity reactions in (/>i?30maj:8- A cycle 
realizes its operation, two cycles reinitialize the beam state and the operations 
can then be repeated. 



Therefore a beam routing to represent such an operation in 0i?3Omaj:8 re- 
quires two particles (as positrons but one slow and other fast respectively) with 
the same orientation and one colhsion contact point on the beam routing. Thus 
Fig. [9] illustrates how a beam routing is designed to produce such periodic col- 
hsions. 




Figure 10: Transition between two beam routing synchronizing multiple reac- 
tions. When the first set of collisions are done a new beam routing is defìned 
with other particles, so that when the second set of collisions is done then one 
returns to the initial condition of the first beam, constructing a meta-glider or 
mesh in (^rhq. 



In this way, we can design more complex constructions synchronizing mul- 
tiple collisions with a diversity of speeds and phases on different particles. Fig- 
ure [10] displays a more sophisticated beam routing design, connecting two of 
beams and then creating a new beam routing diagram where edges represent a 
change of particles and collisions contact point on ECA (^rhq. In such a tran- 
sition, a number of new particles emerge and collide to return to the first beam, 
thus oscillating between two beam routing forever. To better understand this 



doublé beam routing dynamic, consider Fig. 11 where we see multiple collisions 
between particles (first beam routing): 

PaiPa ^ Pb^Pb^Pb 
changing to the set of particles (second beam routing): 

Pm "^Pe^Pè 

defining two beam routing connected by a transition of collisions as: 
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Figure 11: Synchronization of multiple collisions in (pmio on a ring of 1,060 
cells in 1,027 generations, starting with 50 particles from its initial condition. 



(Pa^Pa ^Pb^Pb^Pb) (^^4 ^Pe^Pé)^ and 



^p%,p%) {Pa^Pa ^Pb^Pb^Pb) 

So we see that a beam routing representation not only helps in designing 
collisions but also to implement computation. 

We now explain the function of a cyclic tag system (CTS) in (pmio in beam 
routing terms and discuss a number of theoretical implications. 

It is well known that ECA (pmio is universal [131 |6j. However, a number 
of details about its construction, e.g. a cyclic machine, are uncertain 
Simplifìed implement ations of universal computation in CTS are provided i 

ET]. 



, ci liUlliUei 

;rtain [46 
Àded in ^ 



complete description of CTS working in ^rhq is available at [http : //uncomp . uwe . ac . | 
|uk/genaro/rulellO/ctsRulellO .html| 
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particle 


velocity 


A 


2/3 ^ 0.666666 


B 


-2/4 = -0.5 




-6/12 = -0.5 




-6/12 = -0.5 


Ci 


0/7 = 


C2 


0/7 = 


Cs 


0/7 = 


Di 


2/10 = 0.2 


D2 


2/10 = 0.2 




-4/15 ^ -0.266666 


E 


-8/30 ^ -0.266666 


F 


-4/36 ^ -0.111111 




-14/42 ^ -0.333333 


H 


-18/92 ^ -0.195652 


gun 


-20/77 ^ -0.259740 



Table 2: Properties of particles in (^rhq. 

ECA (^RiiQ has a unique complexity due to great number of particles that 
emerge in the automaton evolution. Table |2] presents a summary of the basic 
particle properties in ìPruq^ following Cook's nomenclature. We can appreciate 
how diverse they are. There are even particles that can expand size forever 
which increases the number of collisions that we can produce on this ECA. 

Initially, a CTS works as a traditional tag system [47 , i.e., as a substitution 
system reading the first symbol on the tape, deleting and putting new symbols. 
CTS are new machines proposed by Cook [13] as a tool to implement computa- 
tions in (fRiiQ. CTS are a variant of tag systems: they have the same action of 
reading a tape in the front and adding characters at the end, nevertheless there 
are some new features as follows: 

1. A CTS needs at least two letters in its alphabet (/i > 1). 

2. Only the first character is deleted (z^ = 1) and its respective sequence is 
added. 

3. If the machine reads a character zero then the production rule is always 
nuli (0 e, where e represents the empty word). 

4. There are k sequences from /i* which are periodically accessed to specify 
the current production rule when a nonzero character is taken by the 
system. Therefore the period of each cycle is determinate by k. 

This way a cycle determines a partial computation over the tape. No par- 
ticular halt conditions are specified in the originai paper [13 , possibly because 
the halting is a direct consequence of tag systems. Let us see some samples of a 
CTS working with /i = 2, /e = 3 and the production rules: 1 ^ 11, 1 ^ 10 and 
1 ^ e. 
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Figure 12: A diagram of a CTS working in ^rhq. 



Our intention is to represent CTS function in ^rhq based beam routing 
of symbols. We employ transition beam routings to explain criticai changes 
that occur due to collisions and sets of particles that are present on each beam 
routing at every instant description of the machine. 

We start with a description of main stages of encoding packages of particles 
in their initial condition. This will be an equivalent of beam routing diagram. 
We will first show how particles and their collisions emulate a CTS in ^Pruq. We 
use packages of particles to represent data and operators. We read, transform 
and delete data in the tape using reactions between the particles. The approach 
is a delicate one and requires a laborious task of setting initial confìgurations 
for the particles in the beams. A diagram of such represent ation is shown in 
Fig. [T2j particular features of the diagram are explain below. 

A construction of the CTS in ipRUQ can be partitioned, essentially, in three 
parts: 



First, is the left periodic part, controlied by packages of A_A^ particles. 
This part is static and controls the production of O's and l's. 

The second part is the center, determining the initial vaine in the tape. 

The third one is the right cyclic part, which has the data to process, adding 
a leader component which specifìes data added to or erased from the tape 
in the evolution space. 



In the left part, the four packages of (Fig. 13 3) particles must be care- 
fully explained because although they are static, their phases change periodi- 
cally. The important point is to implement these components in defìning both 
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distances and phases, because a distinct phase or a distance induces an undesir- 
able reaction; in fact ali components obey this restriction because every glider 
of each component must be correctly aligned. 

The centrai part is constituted by an initial 1 on the tape represented by a 
package of four C2 particles. This way, an element lEle_C2 (Fig.flSb) represents 



a 1 and the element 0Ele_C2 (Fig. 13 1) represents a in the tape. 

The element OBlo_^ is formed by 12È particles. For the lBlo_^ element, 
Cook establishes the existence of two components to represent l's: lBloP_^ 
(Fig. 13:) named primary and lBloS_^ (Fig. 13 5) named standard. In essence 
both blocks produce the same element lAdd_£^, although coming from different 
intervals. The reason to use both blocks is that (pmio is not symmetric; thus, if 
we use only one element then although we would obtain a good production in 
the first collision, generating an element lAdd_^, when the second package of 
4_A^ particles fìnds the second block, the system is completely destroyed. 

A leader element SepInit_£'E' (Fig. 131) is important in order to separate 



packages of data and determine their incorporation into of the tape. 

) and OAdd_^ (Fig. 13 1) are produced by two 



The elements lAdd_£; (Fig. 
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previous different packages of data. An element lAdd_£^ must be generated by 
a block lBloP_^ or by lBloS_^. This way, both elements can produce the same 
element. While an element OAdd_^ is generated by a block OBlo_^ (Fig. 13 3). 
Finally, in terms of periodic phases, this CTS can be written as follow: 

left: {649e-4_A'^(FJ)}*, for 1 < i < 3 in sequential order 

center: 246e-lEle_C2(A,fi_l)-e-A3(fi_l) 

right: {SepInit_£^£^(#,fi_l)-lBloP_£^(#,fi_l)-SepInit_^^(#,fi_l)-lBloP_£^(^ 

OBlo_£;(#,fi_l)-lBloS_^(#,fi_l)}* (where 1 < z < 4 and # represents a 
particular phase). 

For a complete and full description of such reproduction by phases fi_l, 
please see [46j . 

To get a CTS emulation in (pmio by beam routings, we will use connections 



between beam routing, as we have proposed in Fig 10 



Transitions between beam routings connect each set of collisions to enter to 
the next beam routing diagram. 

We need some fine details to get a CTS operation by beam routings. Fig- 
ure 14 shows the general diagram to reproduce a CTS by particle collisions in 



^Ruo with beam routing transitions. While Fig. 13 describes each component 
to code particles in (pmio, Fig. [M] display how to connect such components and 
control the transition of collisions. 

Some notes are necessary to better understand this schematic diagram. For 
the components lEle_C2 and 0Ele_C2 they are compressed only for one dark 
circle (that represents the point of collisions). Both elements are constituted for 
four C2 particles in different distances, although they have a static position 
that can be confìned to this dark circle, as p^ and Pq respectively. 
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Figure 14: Beam routing machine transitions to simulate CTS in (fmio. 
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When a leader component (Sepinit _£'£') is transformed, given previous bi- 
nary value on the tape, it collides with component, i.e. a or p^ element. 
If is 0, then a cascade of collisions start to delete ali component s that come 
with three particles successively and change this reaction to reach a new leader 
component and return to the beginning. But, if p^ is 1 then a cascade of trans- 
formations dominated by additional particles p^ starts, in order to reach the 
next leader component. Here, we have more variants because pre-transformed 
package of particles is encoded to binary values that is then included on the 
machine tape. If a block of particles is lBloP_^ or lBloS_£' this means that 
such a component will be transformed to one lAdd_^ element. But, if a block of 
particles is OBlo_^, then such a component will be transformed to OAdd_^ ele- 
ment. At this stage, when both component s are prepared then a binary value is 
introduced on the tape, a lAdd_^ element stores a 1 (lEle_C2), and a OAdd_^ 
element stores a (0Ele_C2), which eventually is deleted for the next leader 
component and starts a new cycle in the CTS machine. 

The relevant point here is how to control and code particles by beam rout- 
ings. That will offer a better chance to perform a computation with this archi- 
tecture before implementing the complicated and laborious task of setting up 
initial condition. Of course, such an idea can be carried out in any CA handling 
signals, waves, particles, gliders or mobile self-localizations [H [2j [3l [H [121 Eli 

The package particles codifìcation for a CTS in [61] [46] is represented for 
the expression (slsl01)+ while that in JTT is represented for the expression 
(slllsO)^. Finally in ali cases, this CTS beam routing may simulate such 
operations. 

6 Summary 

We advanced a concept of symbol super collider [60] and implemented it in one- 
dimensional cellular automata. Two types of automaton rings were considered 
— elementary cellular automata (binary celi state, ternary neighbourhoods) and 
element ary cellular automata with memory. We demonstrated that, for certain 
rules of cell-state transition functions, the automata exhibit a wide range of 
particles (mobile self-localizations or gliders) in their evolution. High degree 
of morphological diversity and richness of collision dynamics of the particles in 
memory-enriched cellular automata make them particularly promising objects 
for constructing computational schemes. 

In present paper we advanced over twenty years old studies in identifying 
and classifying gliders and localisations in CA. There is a number of approaches 
varying in their degrees of efììciency and discovery potentials. Thus, Wuensche 
successfully classifìed CA evolution rules with his Z parameter, and provided 
power structures of a substantial number of glider-generating rules ^63j. Other 
studies focus on specification of sets of strings as patterns and gliders. For exam- 
ple, Redeker developed a compact algebraic representation of ID CA evolution, 
known as flexible Urne. Suing his approach one can construct specifìc periodic 
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structures [54]. Evolution programming techniques are also proved to be very 
successful, e.g. in situations when a genetic algorithm undertakes a search for 
gliders and patterns on fìxed small arrays of bits [55l|62]. As well, Freire et. al 
developed an approach to select specific sets of vectors, initial configurations, 
that determine development of gliders and localised patterns [18j [19] . Finally, 
the de Bruijn diagrams is a time-tested and reliable tool to classify sets of pe- 
riodic structures or gliders using exhaustive approach [40l [41]. Said that in 
present paper we did not aim to deal with Identification of mobile patterns 
but rather consider glider sets in collision-based computing perspective. Con- 
cretely, we only aimed to develop the idea of cyclotrons as an abstract model of 
super-colliders [60 . 

We proved experimentally that the dynamics of these particles support ali 
basic types of ballistic collisions and consequently that the cellular automata 
cyclotrons can be used to implement certain schemes and operations of collision- 
based computing. We designed schemes with several beams of particles and 
provided a means of beam routing and programming interactions between par- 
ticles. To demonstrate the high degree of computational universality of cellular 
automaton 'cyclotrons' we took into account that each particle is essentially a 
finite-size binary string travelling along a one-dimensional lattice 'cyclotron.' 
When a string with lower velocity overtaken by a string with higher velocity 
the content of one or both strings can be modified. Thus, for carefully selected 
initial configuration of particles and types of particles on a ring, the cellular 
automaton cyclotron imitates a cyclic tag systems, and thus is computationally 
uni ver sai. 

What are advantages of the proposed approach? In cellular-automaton mod- 
els particles can run along the rings indefinitely. With suitable beam routing 
schemes some particles (the results of computation or the products of reactions) 
can be removed from the system and new particles can be added. The cellular- 
automaton super collider is a universal computing devices based on interactions 
between particles due to the particles' different speeds. This give our constructs 
an enormous advantage when compared to more "traditional" constructs, where 
particles collide only due to different orientations of velocity vectors. Also most 
existing unconventional designs of universal computing devices have certain is- 
sues of boundary conditions, e.g. infinite versus finite, the cellular-automaton 
super-colliders by using periodic boundary conditions removes the effect of any 
boundary being equivalent to an infinite periodic system. 

The theoretical models of cellular-automaton super colliders open new per- 
spectives in laboratory implementations of collision based computing devices. 
The opportunities are virtually endless. Hundreds of chemical, physical and 
biological systems exhibit travelling localizations in their dynamics. Examples 
include co-aligned dipole groups in tubulin microtubules (e.g. [531 [39j, kinks, 
breathers and solitons in molecular chains and polymers (e.g. [16], phasons 
in quasi-crystals [28], kinks in ferromagnets [15], dissipative solitons in gas- 
discharge systems [8] , localizations of electron density in monolayers of aromatic 
molecules [10] , wave-fragments in excitable chemical media [4 . We envisage that 
molecular chains and polymers would be the best candidates for experimental 
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laboratory realisation of symbol super colliders. This will be a subject of further 
studies. 
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